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Abstract 

We investigate the effects of fluctuations of the nuclear surface on 
the harmonic oscillator elastic charge form factor of light nuclei, while 
simultaneously approximating the short-range correlations through a 
Jastrow correlation factor. Inclusion of surface-fluctuation effects 
within this description, by truncating the cluster expansion at the 
two-body part, is found to improve somewhat the fit to the elastic 
charge form-factor of ie O and Ca. However, the convergence of the 
cluster expansion is expected to deteriorate. An additional finding 
is that the surface-fluctuation correlations produce a drastic change 
in the asymptotic behavior of the point-proton form factor, which 
now falls off quite slowly (i.e. as const. • q~ 4 ) at large values of the 
momentum transfer q. 
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1 Introduction 



The calculation of the charge form factors F c h(q) of nuclei is a challenging and 
appealing problem A possibility to face this problem is by means of an 
independent-particle model. In this approach, which is particularly attrac- 
tive because of its simplicity, the choice of the single particle potential has 
to be suitably made. In fact a short range repulsion in this potential seems 
advisable for light nuclei (see ref. and references therein). For example, 
with an harmonic oscillator (HO) potential having in addition an infinite soft 
core of the form 4 (B > 0) the F ch (q) of 4 He can be well reproduced, but 
for the heavier nuclei, such as 12 C and 16 0, state dependent potentials seem 
necessary and even then the fit is not so good for higher q- values || . Further- 
more, the correction of the centre of mass motion can not be made exactly 
and unambiguously. An approach which is rather similar is the approach of 
Ripka and Gillespie |§ and Gaudin et al [§]]. It was shown by these authors 
that if a Jastrow wave function consisting of HO orbitals and of simple state- 
independent correlation functions of the form /(r) = (1 — e _/3 r )^ 2 is used, 
one can construct a Slater determinant which yields a density "very similar" 
to that of such a Jastrow wave function in the two-body approximation. This 
is done by diagonalizing the density matrix: 

that is by using the so-called "natural orbitals". It is clear, however, that the 
actual density matrix is unknown and one obtains it approximately, usually 
by means of a Jastrow wave function, as is the case in the above references 
(see expression (7.7) of the paper by Gaudin et al ||]). Therefore, even with 
such an approach one still has in practice to use a complicated wave function 
such as a Jastrow one. From the above discussion it is clear that a Jastrow 
wave function is useful even if one is only interested in calculating the form 
factor and density distribution of nuclei and not necessarily other quantities, 
such as the momentum distribution || . 

In a series of papers || [7], ^] an expression of the elastic charge form 
factor, F c h(q), truncated at the two body term, was derived using the factor 
cluster expansion of Ristig et al || |10|| . This expression, which is a sum 
of one-body and two-body terms, depends on the harmonic oscillator (HO) 
parameter h\ and the correlation parameter A through a Jastrow type cor- 
relation function which introduces the short range correlations (SRC). The 
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use of the HO orbitals (as well as of the particular and popular form of the 
Jastrow correlation function used) stems mainly from the considerable sim- 
plification they imply. There are, however, additional advantages, such as: i. 
The correction of the centre of mass motion can be done exactly by means of 



the Tassie and Barker factor |TT[ . ii. One can obtain analytically the asymp- 
totic behaviour of the form factor F c h(q). In principle, of course, one should 
start with a Hartree-Fock independent-particle model assuming an effective 
interaction and then introduce correlations. Such an approach is, however, 
computation-wise quite demanding. Although it is more sophisticated and 
more satisfactory, it lacks the above mentioned two advantages. 

The fit of F c h(q) to the experimental data with the above mentioned 
procedure was very good both for low and high values of momentum transfer 
except for the values around the last maximum for 16 and m Ca. Better fit 
can be obtained if the parameter A is taken to be state dependent but in this 
case there is a big number of parameters, six for 16 O || and ten for 40 C a 

B 

Another possible way to make the agreement between theory and exper- 
iment better might be to introduce, in addition, other types of ground state 
("long range") correlations which have been the subject of previous investi- 
gations by a number of authors (see, for example, |H| [14], [15], [16|, [H], [18], |19|). 
We focus our attention, as in ref . 1T7TI , on fluctuations of nuclear surface due 



to the zero point motion of collective surface vibrations pO , pif| which can 
affect the ground state charge density. The presence of surface fluctuation 
correlations (SFC) introduce another fitting parameter in addition to the HO 
and the SRC parameters. Thus, it appears to be of interest to develop the 
relevant formalism and to investigate what would be the effect, if any, of this 
additional parameter to the best fit values of the other parameters and to 
the quality of the fit. The aim of this paper is to report on results of certain 
investigations towards this direction. 

In section 2, the above SFC are introduced to the HO densities and an- 
alytic expressions of the nuclear density, elastic form factor and of the n-th 
moment of the density distribution are given for the nuclei 4 He, l6 and 
40 C 'a. In section 3 the introduction of the SFC in addition to the SRC is 
studied. Numerical results are reported and discussed in section 4. 



3 



2 The effect of the surface fluctuations on 
the harmonic oscillator density and form 
factor 

Our starting point is the expression for the proton (or charge) density of 
a nucleus which has been deformed through the zero-point motions of the 
collective surface vibrations. This expression, according to ref. |T7[ (see also 
ref. 



13| , |14| , [15| for a rather similar expression) is the following: 



Pia{r) 



2na 



px(r - f)exp 



2a 2 



(1) 



where pi (r) is the uncorrelated density, s is a correction needed to conserve 
the number of particles in the correlated ground state and a is a measure of 
the effect of the zero point fluctuations. The value of a is related to (3\, the 
deformation parameters for the states of multipolarity A , with the relation 
a 2 ~ "3" X)a Px( t = 0) while the (3\ parameters can be determined from the 
values ofB(E x ) 0, 0. 

In (H) we consider for pi(r) the HO proton density, in which the centre 
of mass correction has been taken into account, for nuclei He to 40 C a. This 
is given by the expression: 
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and b 2 = 6? (1 — 4), A is the mass number and 6i = \/— the harmonic 
oscillator parameter. r\ n i is the occupation probability (0 or 1 in the present 
treatment) of the nl state. It is easily checked that when b\ = b\, that is 
when the centre of mass correction is not taken into account, the coefficients 
of the polynomial in El are reduced to the well known expressions: 



N = 2771s + 3r) 2s , ^2 = 4r] lp - 4r] 2s , N± = -r} U + -r] 2s 
General expressions of similar structure for the density and the form factor 



in the HO model have been given in ref 22 



From (P and (|2|) an analytic expression of pio-(r) can be derived. This is: 
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where the coefficients Ck depend on No, N 2 , N4, a, s and b\ and are given 
by the following formulae: 
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and B = l/(b 2 + 2a 2 ) 

By using expression @ one can find an analytic expression for the nth 
moment of the density. This is the following: 
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An approximate expression for < r n >\ a may be derived by truncation 
of the series at the second power for a and the first power for sq : 
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(7) 



By taking into account that 

< r° >i a =< r° >!= 1 
the approximate expression for the parameter s is: 
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(8) 



That expression was used as a first approximation in our calculations. 
More accurate values were obtained by varying sq until normalization of 
Pio-( r ) was achieved to a good approximation. From expressions (^) and (@) 
and from the known expression of the moments of the HO density one can 
find the approximate expression of the contribution of the SFC, A < r 2 > lcr , 
to the mean square radius for nuclei A He to m Ca. This is given by the 
following expression: 
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Finally, for the elastic point proton form factor the well known expression 
in Born approximation 
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is used. Substitution of pi a (r) from (|J) leads to the following analytic ex- 
pression of Fi a (q) in terms of the confluent hypergeometric function 



where 
h = 
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The complex quantity z is given by: z = vBs + iq/(2yB) . 

Expression ) may be reduced to a somewhat more convenient form: 
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where 



2 C 



1 



A c 2 
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. (14) 

The coefficients C 2n and (72,1 depend also on N , N 2 , a, s and b\ and 
are given by the following expressions: 

C = N + ^N 2 + ^iV 4 + C (2N + N 2 + ^iV 4 ) 

C 2 = -AT 2 - (5 - Aa 2 B)N 4 
C 4 = bjBN 4 

d = § (2iV + N 2 + ^iV 4 ) 

~C 2 = -s 5(2iV 2 + 6iV 4 ) 
^4 = 2s 6 2 5 2 iV4 



(15) 
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It may be easily checked from expression that when the SFC are 
switched off, that is when the limiting case a — > is considered, expression 
( P~5| ) for Fi a (q) goes over to the well known harmonic oscillator one, as should 
be the case, on the basis of expressions (0) and ([[]). Furthermore, by using 
the asymptotic expansion of the confluent hypergeometric function, we find 
that the behaviour for Fi a (q) at large values of the momentum transfer is 
the following: 

F M~ 1 ex P[~ g o B ] ( A ( 1 )-4 , a ( 1 \-6 , a ( 1 \-b\ 

FxM ~^z ijVB v w 6 W 8 W ) 

(16) 

where 

Aa = — so Co — ~2B^ 1 

A, = (-3s + 2Bs 3 )C + (3s 2 -^)C 1 + ^C 2 + ^C 3 

A 8 = + 15S So 3 - 5975 W)C + (~ + ^ + ™Baf)C, + 

r 45s ° i 5s 3 lr , r 45 45s o' r 45 So f r 

{ YW ~ 15 So )C2 + W " Tb" )C3 " 2& Ci (17) 

Thus, it is seen that for sufficiently large values of q, the form factor tends 
to zero rather slowly, namely as the inverse fourth power of the momentum 
transfer. On the contrary, the HO form factor goes rapidly to zero for large q, 
namely as a Gaussian or as a Gaussian times an even power of q (depending 
on the nucleus). 

The value of q at which the Fi a (q) approaches the value given by the 
asymptotic expression ( |T6] ) does not seem to depend very strongly on the 
nucleus, at least when the values of the parameters b\ and a are determined in 
the way described in the following two sections. In Fig. 1 the F\ a (q) has been 
plotted for the 16 O nucleus, using the values &i = 1.563/m and a = 0.414/m 
(see section 4), together with its asymptotic behaviour const. ■ q~ 4 and the 
improved asymptotic expression fliHD, respectively. It is seen that Fi a (q) 
becomes close to the asymptotic behaviour const. ■ q~ 4 at quite large values 
of the momentum transfer ( larger than 10/m _1 ) while essential convergence 
to the behavior fll~6|) is achieved at much smaller q values. It might also be 



of interest to note that people have assumed in the past |23j a decrease 
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of the form factor of the type Cq~ 4 in the region where no measurements 
are performed: q > q max (measured), in order to obtain error envelopes on 
the densities of nuclei. The present analysis indicates that the inclusion of 
additional terms of the type mentioned above seems appropriate in this sort 
of analysis. 

There are two parameters in expression (|I3| ) , the HO parameter b\ and the 
SFC parameter a which can be determined from the deformation parameters 
(3\ associated with the low lying collective states of the nucleus or can be 
treated for example as a free parameter. In the latter case, the fit of the 
form factor ( p!3|) (after correcting it for the finite proton size pi) to the 
experimental data (refs. [24], |25| for A He, |26j for 16 and |27j for 40 Ca ) 



leads to zero value for the parameter a except for A He. For A He, a is different 
from zero (er = 0.706/m and b\ = 1.252/m) and the value of \ 2 ls smaller 
compared to the one obtained with the HO model. However, the diffraction 
minimum is not reproduced. Because of these reasons the introduction of 
short range correlations is advisable. This is done in the next section. 



3 The effect of the surface fluctuations and 
short range correlations on the charge form 
factor and density 

A general expression for the charge form factor F c h(q) of light closed shell 
nuclei was derived || || using the factor cluster expansion of Ristig, Ter Low 
and Clark || 10 1. This formula was further simplified by using normalized 



correlated wave functions of the relative motion which were parameterized 
through a Jastrow type relative two-body wave function of the form: 

Vws(r) = N nls [l - ex P (-\r 2 /b 2 )}<p nl (r) (18) 

where N n i$ are the normalization factors, 0„z(r) the harmonic oscillator wave 
functions and b = \f2b\ is the HO parameter for the relative motion. The 
expression for F(q) is of the form 

F(q) = F 1 (q) + F 2 (q) (19) 
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F\(q) is the contribution of the one-body term to F(q): 



F 1 ( q ) = ^e W [- b -f]^N 2k ( b fr (20) 

k=0 



where 



id) 



N = 2(r] ls + r] 2s + 3r] lp + 5r] ld ) , N 2 = -^(2^ + 3r] lp + 10r] 

(21) 

iV 4 = |(4t7 2s + 877id) 

while F 2 (q) is the contribution of the two-body term to F(q) and is a function 
of q 2 through the matrix elements 

Ani l s S '{jk) =< Tpnis\jk(qr/2)\i) nn , s > > 

It consists of simple polynomials and exponential functions of q 2 . 

The point proton density can be obtained from flI9| ) by Fourier trans- 
forming F(q). The density is also separated into two parts: 

p{r) = pi(r) + p 2 (r) (22) 

Pi(r) and p 2 (r) are the Fourier transforms of Fi(q) and F 2 (q), respectively. 
Pi(r) is given by expression (H) (with bi = b±) while p 2 {r) is calculated 
numerically because of the complexity of F 2 (q) mainly for 40 C a. 

The correlation parameter A and the HO parameter b\ were determined 
by fitting F ch (q) = fp(q)fcM(<l)F(q) to the experimental charge form factor. 

fpio) i fcKiio) — ex P[^"] are ^ ne correc ti° ns due to the finite proton size 
and the centre of mass motion |llj , respectively. 



As it was pointed out in the introduction, a possible way of improving 
the quality of the fit in the approach outlined previously, should be to take 
into account the correlations originating from the fluctuation of the nuclear 
surface. Thus in (f22j), pi(r) is substituted by pio-(r), while p 2 {r) by p 2a (r), 
which results by using in expression ([]]) instead of pi(r) the Fourier transform 
of F 2 (q) = fcM(q)F 2 (q). Therefore, instead of p(r) given by (f22j) we have 
now p a : 

Pa(r) = pia{r) + p 2CT (r) (23) 
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The (total) point-proton form factor F a (q) is then obtained by summing the 
Fourier transforms of pi a (r), F\ a and of p2 (T) F 2a , that is: 

F a = F\ a + F 2a (24) 

It should be noted that because of the assumed form of the correlation 
function the asymptotic behaviour of F a (q) for large q is expected to be of 
the same functional form as that of Fi a (q). The coefficients, however, in the 
various even negative power terms of q will be different. Thus, the introduc- 
tion of the SFC effects in the correlated (through Jastrow correlations), H.O. 
wave function leads to less steep decrease of F a (q) at large q. It might be of 
interest to point out that a slow decrease of the form factor at large q is also 
suggested from our numerical results. 

Finally, the charge form factor is obtained by multiplying F a (q) by the 
charge form factor of the proton, f p (q): 

Fch(q) = f P {q)F a {q) (25) 

Expression ( p5[) was used in fitting to the experimental values of the 
charge form factor, by treating b\, A, and o as fitting parameters. 



4 Numerical results and discussion 

The best fit values of the three parameters in the form factor, as well as the 
values of \ 2 1 f° r the nuclei 16 and m Ca are displayed in table 1 where 
three cases are considered. In case 1 there are no correlations of any kind 
while in case 2 the SRC are included. These two cases have been studied 
in previous works || [7], ||. Finally, in case 3, both the SFC and SRC are 
included. In the case of A He the best fit value of the parameter a goes to 
zero in the least-squares fitting procedure. Therefore, the present approach 
does not seem to be applicable to this nucleus. 

The value a = 0.293 fm for m Ca may be compared with the value a = 
0.638/m which is given in JT8| . It is seen that the value of a is considerably 



smaller in the present approach in comparison to that of ref. JT8|. An 
analogous remark holds for the value of this parameter for 16 0, if comparison 



is made with the value obtained in ref. [30 |. 

It is further noted that the value of a for 16 is larger than that of 40 Ca. 
This is in accordance with what one expects on the basis of the variation 
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of a 2 with the mass number, which was studied in ref. JJIJ and is further 
elaborated in the Appendix of the present paper. We also observe that if 
the parameter k (see expression (A.l)) is determined from our value of a for 



ie O, the value of a we find for 40 Ca is not far from the value predicted on 
the basis of expression ( |A.7|) . It is smaller than that value by less than 30%. 



We observe also from table 1 that for 16 and 40 Ca, the introduction of 
the SFC and the SRC decreases the values of parameters &i and A while the 
value of (^i/A) 1 / 2 (which is the "actual correlation parameter", since small 
values of (fe^/A) 1 / 2 imply values of the correlation factors closer to unity) is 
increased. That would deteriorate the convergence of the cluster expansion. 
This is indicated by the values of the so-called "healing" or "wound" integrals 
(see ref. [[JTJ) for the various states: 



Vnls = / I ^ms{r) - <p n i(r) \ dr (26) 
Jo 

It is expected that, the larger the value of rf the worse the convergence of 
the cluster expansion of the density. With the correlated relative, two-body 
wave function (JIB|), the healing integral is given simply by 



V 2 nls = 2[l + N nls (I nls -l)} (27) 

where the normalization factors N n i s and the integrals I n i s = / °° e~ Xr ' 2 l b2 4> n idr 
can be obtained analytically for the various states. For the lowest state, its 
expression is fairly simple: 



Voo 



2 |l - [l - (1 + A)- 3 / 2 ] [l - 2(1 + A)- 3 / 2 + (1 + 2A)- 3 / 2 ] V2 j (28) 



and can be easily used for calculations. We have dropped the index S, since 
we consider state- in dependent correlation functions. For 16 0, in case 2 we 
find r/ 2 = 0.00717 while in case 3, r/ 2 = 0.01379, indicating deterioration 
of the convergence of the cluster expansion. Similar are the results for the 
higher states (e.g. for the lp state we have in case 2, 77^ = 0.00027 and 
7701 = 0.00083 in case 3). Analogous are the results for 40 Ca. In this case, 
however, the increase of the values of rj^ is smaller in comparison with the 
corresponding one for 16 0, indicating milder deterioration of the convergence 
of the cluster expansion. 

It is clear regarding the results of case 3 in table 1 that the SFC effects 
on both, the one-body and the two-body part of the cluster expansion of the 
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form factor have been taken into account. Quite a simpler approach would 
be to take into account these effects only on the one-body part. Such an 
approach, however, would not be appropriate since then both terms would 
not be treated on the same footing. It might have been acceptable if the 
effect of the SFC on the two-body part had been sufficiently small, which 
unfortunately does not appear to be the case. This is indicated by the change 
in the best-fit values of the parameters. If the SFC effects are taken into 
account only for the one-body part of the form factor, then the best fit 
values for 16 are: 61 = 1.647 fm, A = 11.440, a = 0.224 fm and for 40 Ca: 
b\ = 1.814 fm, A = 11.786, o = 0.364 fm. It is seen that there is a noticable 
change in the values of these parameters. Furthermore, the value of a for 16 
is smaller than the value of o for 40 C a, which contradicts our expectations 
(see Appendix). 

For the two nuclei we have considered, the introduction of the SFC has the 
effect of improving the fit of F c h(q) to the experimental data. Although the 
introduction of SFC does not decrease the value of x 2 too much, its relative 
decrease for 16 being larger than that of 40 C a, there is an improvement of 
the fit in the region of the large g-values. This can be seen in figures 2a 
and 3a where the F ch (q) for 16 O and 40 C a have been plotted with the best 
fit values of the parameters and compared with the experimental F c h(q). 
Another interesting observation one can make regarding these figures is the 
slow decrease of F c h(q) at large q in case 3 (in which the SRC and SFC are 
included in the form factor), in comparison with cases 1 and 2. 

It seems also appropriate to point out that the putative roles of the mean- 
field, short range correlations and surface-fluctuation effects get mixed up to 
some degree. Thus, (see table 1) the introduction of SRC changes the value 
of the HO parameter and the additional introduction of SFC not only does 
it change further this value but also the value of the parameter determing 
the short range correlations. Furthermore, the SFC effects which should 
predominantly influence the low and medium q behaviour of the form factor 
produce a substantial change in its values and its functional behavior at very 
large values of the momentum transfer. One could perhaps say that the 
introduction of the SFC "relieves the burden" assumed by the SRC (in the 
two-body approximation) in correcting the independent particle model, so 
that they can better perform the function for which they were designed. 

In figures 2b and 3b the (corresponding to F c h(q)) charge densities p c h( r ) 
of 16 and 40 C a have been plotted. In the same figures certain other relevant 
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quantities have also been plotted, namely: 

a. The two-body part of the density p 2 (r) corrected for the finite proton 
size: p' 2 (r) without SFC. The dashed line p' 2 (r) was obtained with the pa- 
rameters of case 2 and the solid line p 2 (r) was obtained with the parameters 
of case 3. It is seen that if the two-body part of the cluster expansion of 
the charge density with the Jastrow correlation function and HO orbitals is 
calculated with the parameters of case 3, its absolute values for various r 
increase rather considerably (see solid line) in comparison with its absolute 
values calculated with the parameters of case 2. 

(3. The difference of the charge densities Ap c h(r) = p c h(HO + SRC + 
SFC) - p ch (HO + SRC), that is the difference between the (HO+SRC) 
charge densities with and without SFC, calculated with the parameters of 
case 3. It is seen that Ap ch (r) is quite small and characterized by oscillations. 

In summary, the present analysis suggests that the inclusion of the cor- 
relations originating from the fluctuations of the nuclear surface in the usual 
cluster expansion (truncated at the two-body part) of the charge form fac- 
tor of w O and m Ca leads to some improvement in the quality of the fit to 
the experimental data. It is also found that the convergence of the cluster 
expansion is expected to deteriorate. Furthermore, the inclusion of these cor- 
relations has a drastic effect on the asymptotic behaviour of the point-proton 
form factor, which now falls off for large q quite slowly, that is as const. q~ 4 . 
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A Appendix 



In this appendix we discuss the variation of the parameter a 2 with the mass 
number A on the basis of a phenomenological analysis. We assume that cr 2 is 
proportional to the ratio: (surface of the nucleus) / (volume of the nucleus). 
This is suggested by the results and discussion of ref. pQ| . Thus, we may 
approximately write: 



a 



R 



(A.l) 



eq 



where k is a proportionality constant and R eq the equivalent uniform radius 
of the nucleus, defined by 

3 



< r > eq = -R; q =< r > 



(A.2) 



We consider that the nuclear density distribution is approximated by a 
symmetrized Fermi function |32j: 



PsF(r) 



exp 



r — c 



exp 



— r — c 



(A.3) 



where c is determined by the normalization condition and given by [[33], [34[] : 
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243 Vr A^ 
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(A.4) 



The parameter r is expressed in terms of po : r o = (i^) 
Thus, we may write: 



^2^ ^2^ °2 

< r >~< r >sf= ~c 
5 



(A.5) 



On the basis of expressions (|A.1|) , (|A.2|) and ( |A.5| ), we find: 



a 



1 + 



7raV 
c J 



-1/2 



(A.6) 
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Except for the very light nuclei, we may use only the first terms of the 
expansion of the above expression in powers of A and we may therefore write: 



a 2 ~ kr^A" 1 ' 3 



5 / na 
1 ~ 6 Vr AV3 



71 / na \ 4 

72 \r A^J 



(A.7) 



The parameters r and a appearing in c may be determined by a least 
squares fit of the RMS radius of the symmetrized Fermi distribution to the 
experimental values of the RMS radius of nuclei: < r 2 > X J 2 V - Considering the 
nuclei of table 2 of ref. |3(J and the corresponding values of < r 2 > X J 2 V cited 
there, we find r = 1.147 fm and a = 0.507 fm. The quality of the fit is very 
satisfactory. 

The value of k may also be determined by fitting to the values of a 2 

go 



shown in Table 2. The corresponding best fit value is k 



obtained in 
2.028 fm 3 . 

Our results for the RMS radii are displayed in table 2, where also the 
results for a 2 obtained with (|A.6|) and ( |A.7|) are given, as well as those with 
the expression a 2 = kiA~P. The best fit values of the parameters in the 
latter expression are k\ = 1.415 fm 3 and j3 



= 0.312, which are rather close to 
1.64t4 -0 ' 35 ), shown also in table 



the corresponding values of ref. 
2. 

We finally observe that expression ( |A.7|) gives results very close to those 
obtained with ( |A.6|) . There is only a small difference for the lighter elements. 
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Table 1: The values of the parameters b\, A, (bf/X) 1 ^ 2 and a, the x 2 ar| d the 
RMS radius and the contributions to it from the HO density and the various 
correlations for nuclei 16 and 40 Ca found by fitting to the experimental form 
factors. Case 1 is referred to the HO form factor, Case 2 when the SRC are 
included in the form factor and Case 3 when the SRC and SFC are included in 
the form factor. All the quantities except A and \ 2 are i n f m - 
a: pL b: pTD 



Case 


Nucleus 


bi 


X 


V x 


17 


x 2 


total 


< r 2 , 

ch 

HO 


>l/2 

SRC 


SFC 


Exper. 


1 


16 Q 


1.786 








9013 


2.728 


2.728 








2 


16 Q 


1.679 


12.768 


0.470 




6226 


2.659 


2.577 


0.654 






3 


Wq 


1.563 


7.989 


0.553 


0.414 


6002 


2.655 


2.433 


0.838 


0.655 


2.728 a 


1 


40 Ca 


1.950 








26847 


3.439 


3.439 








2 


i0 Ca 


1.860 


13.915 


0.499 




19930 


3.420 


3.289 


0.936 






3 


i0 Ca 


1.849 


10.356 


0.575 


0.293 


19588 


3.505 


3.272 


1.165 


0.484 


3.482 b 



Table 2: The values of RMS radii (in fm) and of the parameter a 2 (in fm 2 ) 
for a number of nuclei using the results of ref.[30] and the phenomenological 
expressions of the Appendix 



A 


<r 2 > V2 


< r 2 > 1/2 


cr 2 Ref. H 


<x 2 Ref. Eol 


a 2 














(Table 1) 


(1.64A - * 3 *) 


(1.415A- 0312 ) 


expr. 


A. 6 


expr 


A. 7 




16 


2.73 


2.76 


0.576 


0.621 


0.597 


0.570 


0.588 


40 


3.49 


3.43 


0.487 


0.451 


0.449 


0.458 


0.460 


90 


4.26 


4.29 


0.360 


0.340 


0.348 


0.366 


0.367 


208 


5.50 


5.50 


0.235 


0.253 


0.268 


0.286 


0.286 
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Figure Caption 

Fig. 1. The elastic point form factor in the HO model with SFC: F\ a (q) 
for 16 with 61 = 1.563/m and a = 0.414/m (solid line) and its asymp- 
totic behaviour const. q~ A (long dashed line) together with the values of the 
asymptotic expression (|16|) (short dashed line). 

Fig. 2. The charge form factor (3a) and density distribution (3b) of 16 0. 
The experimental points of the form factor are from ref. [26|, while for the 
density from refs. jSS), ^9|. (For various cases and notation see text ). 



Fig. 3. The charge form factor (4a) and density distribution (4b) of 
40 Ca. The experimental points of the form factor are from ref. p?fl , while 
for the density from refs. p8j , ^9 . (For various cases and notation see text ). 
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